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I • Development  of  Optical  Models 
— 

^Optical  modeling  of  hollow  glass  sphere  aerosols  was  undertaken  with  a 
two  stage  approach.  The  first  stage  was  to  develop  a model  for  homogeneous 
or  solid  spheres.  This  was  a logical  means  of  proceeding,  since  modeling 
of  homogeneous  spheres  is  by  far  the  easier  of  the  two  tasks  and  is  the 
area  which  has  received  by  far  the  greatest  attention  from  other  investigators, 
hence  providing  copious  examples  with  which  to  check  our  results.  The  homo- 
geneous glass  sphere  model,  based  on  the  formulation  of  Mie  (1908),  was 
developed  without  any  serious  problems  and  was  successfully  tested, against 
examples  contained  in  Wickramasinghe  (1973) , Deirmendjian  (1969)  and  Kuriyan 
et  al.  (1974). 

Figures  1-5  depict  a sample  of  the  results  obtained  from  the  optical 
model  for  homogeneous  spheres.  The  quantities  which  appear  in  Figures  1-5 


and  those  pertinent  to  the  optical  model  are 


vext 


Cext/,'lr- 


^ext  *^scat  + ^abs’ 

"ext  = Oext 
x = 2ttR/A, 


(1) 

(2) 

(3) 


(4)  , 

m = n - i k , (5) 

where  Cext  is  the  extinction  cross-section  of  a single  aerosol,  R is  the  aerosol 
radius  and  Qgxt>  Q$ca>  QabS  are  ^e  dimensionless  extinction,  scattering 
and  absorption  efficiencies  respectively.  The  dimensionless  size  parameter 
is  denoted  by  x,  the  wavelength  of  illuminating  radiation  by  A and  the  complex 
refractive  index  of  the  aerosol  by  m.  Figure  1,  using  m = 1.51  - O.Oi,  is 
representative  of  homogeneous  spheres  made  of  ordinary  crown  glass.  The 
range  of  aerosol  sizes  used  in  this  calculation  corresponds  to  the  measured 
range  of  sizes  for  the  Emerson  and  Cuming  FTF-15  aerosols  (henceforth referred 
to  as  FTF-15  aerosols)-  Figure  2,  with  m = 1.5.  - O.li,  depicts  the  behavior 
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of  a similar  glass  with  an  absorbing  dye  added.  Figures  3-5  are  included 
to  show  the  increase  in  spacing  between  extinction  peaks  which  result  from 
decreasing  the  refractive  index. 

^The  second  stage  of  development  of  optical  models  consisted  of  the 
much  more  difficult  task  of  developing  a model  for  composite  or  hollow 
spheres^The  original  formulation  of  the  solution  as  given  by  GCttler 
(1952),  also  appearing  in  k'ickramasinghe  (1973),  was  found  to  be  accurate 
and  stable  only  for  aerosols  whose  product  of  radius  and  refractive  index 
fell  below  a certain  value.  For  aerosols  sizes  and  refractive  indices  in 
the  range  of  interest  for  the  problem  under  consideration,  the  equations 
were  found  to  be  numerically  unstable.  The  equations  were  rewritten  in  a 
numerically  stable  but  analytically  equivalent  form  and  new  give  satisfactory 
results  for  all  ranges  of  input  parameters.  The  model  was  tested  against 
the  published  results  of  Battan  et  al.  (1970),  Wickramasinghe  (1973)  and, 
in  the  limit  of  shell  and  core  materials  of  the  same  refractive  index,  against 
the  results  of  the  homogeneous  sphere  model. 

jgures  6-19  represent  a sample  of  results  from  the  optical  node!  for 
composite  spheres.  The  modeled  aerosol  consists  of  two  concentnc  spheres 
composed  of  two  different  material^,  Notation  is  consistent  witn  that  of  the 
'earlier  model  with  the  addition  of  m^  and  m?,  which  denote  the  refractive 
indicies  of  the  core  (inner  sphere)  and  shell  (outer  sphere),  respectively. 

The  core  radius  is  denoted  by  RQ  and  the  shell's  outer  radius  oy  R,  where 
R-Rq  is  the  shell  thickness  in  pm.  Figures  5-19  show  the  ne  avior  of  a 


hollow  glass  sphere  aerosol,  composed  of  fused  quartz.  The«e  results  are 

'O' 

very  different  in  appearence  from  those  of  the  solid  glass  spheres.  Qualita- 
tively, this  departure  can  be  understood  as  a combination  of  two  phenomena, 
a solid  sphere  of  variable  refractive  index  at  large  aerosol  radii  and  a 

thin  film  interference  phenomenon  at  small  aerosol  radii.  The  long, 

\ 
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slowly-damped  tails  and  lack  of  extinction  peaks  for  large  aerosol  radii  can 
be  simulated  by  a solid  glass  sphere  whose  index  of  refraction  decreases  with 
increasing  size.  Figures  20  and  21  show  the  results  of  such  calculations, 
where  the  indices  of  refraction  were  taken  to  be 

m = (1.45)(R-R0)  + I.OO(Rq)  +0.0i,  (6) 

R 

with  R-RQ  = 1.0pm,  2.0pm  for  Figures  20  and  21  respectively.  Equation  (6) 
merely  represents  a linearly  averaged  refractive  index  for  a hollow  glass 
sphere.  The  behavior  exhibited  in  Figures  20  and  21  can  easily  be  understood 
as  an  extention  of  the  results  presented  in  Figures  3-5.  Tne  resonant  minima 
and  maxima  in  Qgx^.  (Figures  6-19),  for  small  aerosol  radii,  can  be  thought 
of  as  a thin  film  interference  phenomena.  The  observed  dependence  on  shell 
thickness  leads  to  this  hypothesis,  however  the  dependence  can  by  no  means 
be  completely  described  by  simple  thin  film  interference  formula. 
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II.  Emerson  and  Cuming  aerosols,  physical  measurements 

Examination  of  the  physical  properties  of  the  FTF-15  aerosols  began 
with  a preliminary  investigation  using  the  University  of  Wyoming's  Scanning 
Electron  Microscope  facility.  This  investigation  was  intended  to  determine 
whether  their  size  distribution  and  shell  thickness  could  accurately  be 
measured  with  such  an  instrument  or  whether  different  methods  needed  to  be 
explored.  A decision  was  reached  that  the  scanning  electron  microscope 
was  well  suited  for  accurate  determination  of  size  distribution  and  for 
approximate  determination  of  shell  thickness. 

Figure  22  is  a photograph  taken  through  the  electron  microscope,  of 
FTF-15  aerosols  mounted  on  a suitable  substate.  This  photo  is  one  of  a 
series  taken  expressly  for  size  determination  purposes.  A magnification 
of  300  was  used  on  the  electron  microscope,  followed  by  additional  photo- 
graphic enlargement;  a 15pm  scale  is  indicated  on  the  Figure.  Mechanical 
counting  and  sizing  of  such  photographs  lead  to  the  size  distribution  shown 
in  Figure  23.  Shown  is  the  raw  data,  taken  with  lym  resolution  in  diameter 
Figures  24  and  25  are  the  same  data  smoothed  to  2ym  and  5um  resolution 
respectively.  The  mode  of  the  distribution  occurred  at  13ym  diameter,  the 
mean  at  18ym  and  the  spread  from  5 lo  r aim.  The  distribution  rises  sharply 
from  the  minimum  to  the  mode  diameter  then  exhibits  a slow,  exponential -like 
tail  at  large  sizes. 

Figure  26  is  an  electron  microscope  photograph  of  FTF-15  aerosols  after 
they  have  been  crushed  with  a mortar  and  pestal . The  microscope  magnifica- 
tion was  1000,  followed  by  additional  photographic  enlargement;  a 4ym  scale 
is  indicated  on  the  Figure.  With  the  aid  of  this  and  similar  photographs, 
an  upper  limit  of  lym  could  be  placed  on  the  aerosols'  shell  thickness.  The 
aerosol  fragments  are  all  fragments  of  spherical  shells.  Hence,  even  if  viewed 
edgewise,  their  apparent  thickness  will  always  exceed  the  actual  shell  thick- 


ness.  Consequently,  this  technique  can  ,nly  give  an  upper  limit  to  shell 
thickness. 
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A second  and  more  accurate  means  of  determining  shell  thickness,  involvim 

density  measurements,  was  devised.  First,  a few  cm  sample  of  aerosols  was 

suspended  in  water,  in  a hypodermic  syringe,  and  the  combined  volume  measured. 

A hypodermic  syringe  was  needed  because  unless  the  aerosol  liquid  suspension 

was  put  under  some  pressure,  the  aerosols  would  float  atop  the  liquid 

rather  than  be  submerged  in  it.  The  volume  of  the  water  and  the  mass  of  the 

aerosols  were  then  determined,  allowing  for  calculation  of  the  true  average 

3 

density  of  the  aerosols.  Multiple  trials  yielded  .268  g/cm  ± 2%  for  the 
FTF-15  hollow  glass  spheres.  Knowing  the  true  density  of  the  bulk  material 
of  which  they  are  composed,  allows  for  calculation  of  a mean  wall  thickness. 

By  this  method,  the  mean  wall  thickness  of  the  FTF-15  hollow  glass  spheres 
is  0.45ym,  not.  at  variance  with  the  upper  limit  of  lym  from  microscopic 
methods. 

The  final  input  parameter  required  for  the  optical  model  is  that  of 
refractive  index.  The  manufacturer  was  unable  to  give  us  a bulk  sample  of 
glass  for  measurement,  as  the  glass  is  mixed  in  the  liquid  state  and  then 
directly  formed  into  aerosols.  The  manufacturer  however,  seemed  somewhat 
certain  that  the  glass  was  similar  to  boro  silicate  crown,  which  has  a re- 
fractive index  of  1.51-J.0i  at  1.06ym.  Consequently,  we  used  this  value 
in  all  of  our  early  work.  Later,  however  we  received  a private  communi- 
cation from  Huffman  (University  of  Arizona)  stating  quite  conclusively, 
based  on  melting  tests  and  refractive  index  measurements  in  the  visible, 
that  the  FTF-15  aerosols  were  composed  of  fused  quartz.  Consequently,  we 
used  the  index  for  fused  quartz,  1.45-0.0i  at  1.06ym,  for  the  final  calcu- 
lations included  in  this  report. 
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III.  Applications  of  the  Optical  Model 

Combining  the  previously  measured  size  distribution,  the  ,45pm  shell 
thickness  and  the  refractive  index  of  1 . 45-0. Oi » we  are  now  able  to  calcu- 
late volume  extinction  coefficients  for  the  FTF-15  aerosols.  The  volume 
extinction  coefficient,  £ is  defined  as 


8ext  = N-10"3-  t rr2  Qext(m,2Trr/A)C  n(r)dr,  (7) 

-1  -3 

where  Bgxt  is  in  km  , N the  aerosol  number  density  in  cm  , r the  aerosol 

radius  in  pm  . Attenuation  of  a collimated  light  source,  such  as  a laser, 

can  then  be  calculated  simply  from  Beer's  Law, 


F = F0  exp{-0ext  - L>  (8) 

p 

•.,'here  F0  is  the  flux  density  (W/cm  ) incident  upon  the  medium,  F the  flux 
density  exiting  the  medium  and  L the  length  of  traverse  in  the  medium.  If 
the  optical  properties  of  the  medium  are  not  homogeneous,  then  the  simple 
Bext  • L product  is  replaced  by  an  integration  of  Bext  over  the  path  traversed. 

Table  1 gives  the  volume  extinction  coefficient  for  the  FTF-15  aerosols, 
for  various  values  of  suspended  number  densities  N. 

Table  1 

N(cm~3)  0ext^km"^  N (cm-3 ) 3ext(km"‘) 


.01 

1.22  x 

loi 

5.00 

6.09 

.05 

6.09  x 

10- 
10  1 

10.00 

1.22 

X 

10 

.10 

1.22  x 

50.00 

6.09 

X 

10, 

.50 

1.00 

6.09  x 
1.22 

100.00 

1.22 

X 

10' 

Bext  of  course  just  scales  as  a multiple  of  N,  since  N is  just  a constant 
which  multiplies  the  integral  in  Eg.  (7).  To  obtain  a feeling  for  aerosol 
number  densities,  one  should  note  that  naturally  occurring  stratospheric 

_3 

aerosols  have  typical  number  densities  of  1.0cm  and  naturally  occurring 

o _ o 

tropospheric  aerosols  have  number  densities  of  t 10  cm  . However,  these 
naturally  occurring  aerosols  have  volume  extinction  coefficients  of  only 
10‘3  - 10”4  km’1  and  10’ ' - 10"3  km"^  respectively.  The  descrepancy  in 
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the  ratio  of  8ext/N  for  naturally  occurring  versus  FTF-15  aerosols  stems  from 
the  fact  that  the  FTF-15  aerosols  are  orders  of  magnitude  larger  in  size,  the 
factor  r being  very  prominent  in  Eg.  (7). 

Certain  practical  questions  regarding  the  implementation  of  the  FTF-15 

can  easily  be  answered  from  the  results  of  Table  1.  For  instance,  if  it  were 

-1  3 15 

desired  to  produce  a ggxt  of  - 1.0km  over  a volume  of  1km  , a total  of  10 

3 

aerosols  would  have  to  be  approximately  evenly  dispersed  over  the  1 km  volume. 

The  second  primary  application  of  the  optical  model  was  to  investigate  the 
question,  whether  other  types  of  aerosols  might  be  a better  or  somehow  optional 
choice  for  producing  extinction  of  1.06ym  radiation.  To  this  end,  the  optical 
model  was  used  to  calculate  extinction  coefficients  for  aerosols  of  differing 
size  distributions,  shell  thickness  and  refractive  indices. 

Regarding  the  question  of  varying  size  distributions,  one  could  make  a de- 
tailed, but  unrealistic  study  of  various  monodisperse  aerosols.  However,  truly 
monodisperse  aerosols  are  difficult  or  impossible  to  produce  by  high  volume  manu- 
facturing processes.  Hence,  I choose  to  take  the  measured  FTF-15  aerosol  size 
distribution  and  shift  it  to  larger  or  smaller  sizes  by  multiplying  each  size 
bin  by  some  constant  S.  This  procedure  produced  a distribution  which  maintained 
the  original,  realistic  shape,  but  contained  larger  or  smaller  mean  diameters. 
Specifically,  if  the  original  distribution  had  a mean  diameter  of  18ym,  a 
spread  of  5-51  pm  and  the  constant  S were  given  a value  of  0.1,  then  the  new  dis- 
tribution would  have  a mean  diameter  of  1.8ym  and  a spread  of  0.5-5. lym.  Table  2 

_3 

iresents  the  results  of  these  calculations,  where  N = 1cm  was  used. 

Table  2 


s eext(km"1) 

.01  3.35  x 10"? 
.03  5.55  x 10": 
.10  1.48  x 10"f 
.30  1.09  x 10 


5 6ext(km-’) 

1.00  1.22  , 

3.00  1.13  x 10, 

10.00  1.27  x 10^ 


The  results  are  as  to  be  expected,  larger  aerosols  yield  substan- 
tially greater  extinction,  for  a fixed  number  density.  Furthermore, for 


large  sizes  such  that  most  of  the  aerosols  are  to  the  right  of  the  first 

extinction  peak  (i.e.,  see  Figures  1-21),  8gxt  is  approximately  propor- 

2 

tional  to  the  mean  radius  squared  which  is  directly  proportional  to  S . 
However,  larger  aerosols  weigh  more  and  constitute  a larger  volume  of 
agent  to  deliver  to  the  designated  area.  Hence,  if  the  capacity  of  the 
delivery  system  is  included,  there  must  be  some  tradeoffs,  to  consider. 
Let  us  define  a quantity  p. 


p = N ■ 10 


-12 


r°° 


i11  [r3  - (r  - .45pm)3]  • pM  • cn(r)dr,  (9) 
o J 9 


, 3 

where  p is  the  density  (g/cm  ) of  the  aerosol  in  its  suspended  state,  pg 

_3 

the  density  of  fused  quartz,  and  N = 1cm  . Hence  for  a given  volume  into 
which  the  aerosol  is  to  be  dispersed,  p is  proportional  to  the  mass  of 
aerosol  which  the  delivery  system  must  be  capable  of  handling.  Table  3 
presents  the  ratio  of  Bext/p  versus  S. 


Table  3 


S Bg^km'Vpte/cm3)  S Bext(km_1  )/p(g/cm3) 


.01 

2.59  x 

10S 

1.00 

9.24  x 

io: 

.03 

1.60  x 

io| 

3.00 

9.29  x 

io: 

.10 

1.63  x 

,09 

10.00 

9.32  x 

10' 

.30 

1.01  x 

109 

The  units  and  absolute  magnitudes  here  are  not  of  great  importance, 
but  the  relative  magnitudes  are.  These  calculations  show  that  an  S of 
around  .10  would  yield  the  maximum  optical  extinction  for  a given  mass 
of  aerosol  delivered  to  the  designated  area.  An  S of  .10  corresponds  to 
a size  distribution  with  a mean  diameter  of  1.8pm  and  a spread  of  .5  to 
5.1pm.  It  should  be  noted  that  the  small  aerosols  produced  by  S values 
of  less  than  unity  have  been  assumed  to  be  solid,  rather  than  hollow, 
when  their  radius  decreased  to  less  than  .45pm.  This  assumption  is  con- 
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sistent  with  a constant  shell  thickness  under  conditions  of  a variable 
size  distribution. 

The  results  of  Table  3 can  be  easily  understood  if  one  makes  use  of 

2 

the  fact  that  p is  approximately  proportional  to  r for  large  S values, 

3 

and  to  r for  small  S values.  Also  6gxt  is  approximately  proportional 
to  r for  large  S values  and  proportional  to  r for  small  S values.  This  ex- 
plains the  observed  behavior  that  0ext/P  is  nearly  constant  for  .3  < S < 10.0, 
and  drops  off  to  its  smallest  values  at  S = .01. 

A second  means  of  including  delivery  system  capacity  in  tradeoff  de- 
cisions, is  to  calculate  the  packed  volume  of  the  aerosols  which  need  to  be 
delivered  to  the  designated  area.  Since  the  packed  aerosols  are  of  such 
small  density  (.16g/cm  ),  the  volume  rather  than  mass  might  be  a more  im- 
portant quantity  to  consider. 

To  this  end,  we  define  a quantity  V, 


V = • N • 10-12. 


4 tt  r3  • C n(r)dr, 
o J 


(10) 


which  is  dimensionless  and  equal  to  the  packed  volume  of  aerosol  per  unit 

_3 

of  suspended  volume.  The  number  density  N remains  1 cm  and  the  ratio 
.268/. 16  represents  the  ratio  of  true  (liquid  immersion)  to  packed  den- 
sity of  the  FTF-15  aerosols.  Hence,  for  a fixed  volume  over  which  the 
aerosol  is  to  be  dispersed,  the  ratio  of  3ext/V  gives  the  optical  extinc- 
tion per  unit  volume  of  delivery  system. 


Table  4 


s 

BextO^/V 

S 

W®  ); 

.01 

2.71  x 1 0q 

1.00 

9.87  x 10- 

.03 

1.67  x 10q 

1.20  x 10* 

3.27  x 10° 

3.00 

3.40  x 10: 

.10 

.30 

10.00 

1.03  x 10 

These  results  show  a greater  variance  with  S than  for  Table  3,  re- 

3 

suiting  from  the  fact  that  V is  directly  proportional  to  r for  all  values 
of  S.  The  results  show  that  S values  of  .03  to  .10,  a substantially 
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smaller  aerosol  than  the  FTF-15,  would  be  the  optional  choice  if  the  volume 
which  could  be  carried  by  the  delivery  system  was  a limiting  factor. 

A second  major  consideration  in  searching  for  an  optional  aerosol 
for  the  extinction  of  1.06pm  radiation  is  that  of  varying  the  shell  thick- 
ness of  the  hollow  glass  spheres.  Table  5 presents  the  results  of  this 

_3 

study,  where  N was  again  fixed  at  1.0cm  . 

Table  5 


R-Ro(pm) 

6ext(km'1 1 

R-Ro(pm) 

6ext(kl 

1.25 

.771 

.53 

.992 

1.06 

.378 

.50 

1.08 

1.00 

.464 

.45 

1.22 

.79 

.650 

.40 

1.19 

.75 

.791 

.35 

1.06 

.60 

.957 

.26 

.834 

.55 

.959 

.25 

.802 

These  results  show  that  our  measured  mean  shell  thickness  of  .45pm,  for  the 
FTF-15  aerosols, is  approximately  tuned  to  the  1.06pm  wavelength  and  maxi- 
mizes extinction.  The  results  presented  in  Table  5 follow  directly  from 
the  Qext  calculations  shown  in  Figures  6-19.  It  should  be  pointed  out 
that  variations  in  shell  thickness  produce  a variation  of  at  most  a factor 


of  three  in  6ext»  whereas  variations  in  size  distribution  produced  orders 
of  magnitude  changes  in  extinction. 

A third  major  consideration  in  optimizing  aerosol  extinction  is 
that  of  producing  FTF-15  like  aerosols  composed  of  a different  substance. 
To  that  end,  Table  6 shows  the  effect  of  varying  the  refractive  index. 


m. 


Table  6 


1.00- 0. 0i 

1 .00- 0. Oi 

1 .00- 0.01 

1 . 00-0. Oi 

1 . 00-0. Oi 

1 . 00-0. Oi 

1 . 00- 0. Oi 

1.00- O.Oi 

1 . 00- 0. Oi 


1 . 33- 0. Oi 
1 . 40-0. Oi 
1 . 45-0. Oi 
1 . 50-0. Oi 
1 . 55-0. Oi 
1.60-0.0i 
1 . 65-0. Oi 
1 . 70-0. Oi 

1.33- O.li 


ext'1®'1) 

1.07 

1.20 

1.22 

1.15 

1.02 

.89 

.83 

.78 

.90 
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Table-  6 continued 


ml 

m2 

6ext(km 

l.OO-O.Oi 

1 . 70-0. 1 i 

.77 

1.33-0.0i 

1 . 45-0. Oi 

.73 

1-45-0- Oi 

1.45-C.0i 

.74 

1 . 70-0. Oi 

1 . 45-0. Oi 

.74 

1.33-0. li 

1 . 45-0. Oi 

.73 

_3 

Number  density  again  was  fixed  at  1.0cm  . A refractive  index  of  1.33- 
O.Oi  corresponds  to  water,  and  m = 1.70-0.0i  corresponds  to  very  dense 
flint  glasses.  Imaginary  indices  greater  than  zero  correspond  to  adding 
absorbing  dyes  to  the  basic  substance.  Surprisingly,  the  refractive  index 
of  fused  quartz  seems  approximately  tuned  to  the  1.06ym  radiation  just 
as  was  the  .45ym  shell  thickness.  However,  varying  the  refractive  index 
of  both  core  and  mantle  has  an  even  smaller  effect  on  changing  the  extinc- 
tion than  did  varying  the  shell  thickness.  Hence  aerosol  composition  is 
the  least  important  of  the  variables  studied. 

In  concluding  the  study  to  optimize  aerosol  extinction  at  1.06ym,  we 
must  remark  that  we  have  by  no  means  conducted  an  exhaustive  study.  Clearly 
there  are  endless  combinations  of  variables  which  were  not  tried.  However, 
we  did  allow  each  of  the  input  parameters  to  the  optical  model  to  vary  in- 
dependently of  each  other  and  studied  the  response  regarding  extinction. 

We  conclude  that  size  distribution  variations  have  the  largest  impact  on 
extinction,  with  shell  thickness  and  refractive  index  placing  a distant 
second  and  third  respectively.  For  a given  number  density,  extinction 
increases  proportional  to  the  average  size  of  the  aerosol.  However, 
if  one  is  concerned  with  either  the  weight  or  volume  capacity  of  the  delivery 
system,  then  aerosol  distributions  with  mean  diameters  of  0.5-1.8ym  (10-30 
times  smaller  than  the  FTF-15)  are  the  optimum  choice.  The  optimum  size 
is  much  more  sharply  defined  for  volume  considerations  than  for  weight  con- 
siderations. 
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